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ABSTRACT 

We report on Chandra X-ray and Parkes radio observations of IGR Jl 1014-6103, which is a possible pulsar 
wind nebula with a complex X-ray morphology and a likely radio counterpart. With the superb angular resolu- 
tion of Chandra, we find evidence that a portion of the extended emission may be related to a bow shock due to 
the putative pulsar moving through the interstellar medium. The inferred direction of motion is consistent with 
IGR Jl 1014-6103 having been born in the event that produced the supernova remnant (SNR) MSH 1 1-61 A. 
If this association is correct, then previous constraints on the expansion of MSH 11-61 A imply a transverse 
velocity for IGR J11014-6103 of 2,400-2,900 km s"', depending on the SNR model used. This would surpass 
the kick velocities of any known pulsars and rival or surpass the velocities of any compact objects that are 
associated with SNRs. While it is important to confirm the nature of the source, our radio pulsation search did 
not yield a detection. 

Subject headings: stars: neutron — X-rays: stars — pulsars: general — ISM: supernova remnants — stars: 
individual (IGR Jl 1014-6103, MSH 1 1-61 A, SNR G290. 1-00.8) 



1. INTRODUCTION 

IGR J 11014-6103 was discovered as a hard X-ray (20- 
lOOkeV) source tha t was first seen during INTEGRAL 
( Winkler et al.l l2003h observations of the Galactic Plane 
Birdetal.1120101) . Approximatelv 500 "IGR" source^ have 
been discovered or detected for the first time in hard X-rays 
by INTEGRAL. In the Galactic Plane, the INTEGRAL posi- 
tions, which are typically accurate to a few arcminutes, are 
usually not sufficient to obtain source identifications. Thus, 
the field was followed-up with Swift observations, and a soft 
X-ray counterpart was detected (Malizi a et al.ll20ril) . 

Pavanetal. (2011) report on ROSAT, ASCA, Einstein, 
Swift, and XMM-Newton observations of the region. While 
IGR Jl 1014-6103 was the target for Swift, the other satellites 
observed the region in order to study the supernova remnant 
MSH 1 1-61A (SNR G290 .1-00.8). From the X-ray obser- 
vations, [P^inetin daonb find that IGR Jl 1014-6103 has a 
complex X-ray morphology, consisting of a point source, an 
extended source that is ~'22" North-East of the point source, 
and a streak of emission that exten ds to the North- W est of 
the point source and is ^4' in length. iPavan et aTl (1201 1 ) also 
search optical/IR, radio, and gamma-ray catalogs for coun- 
terparts, and they report a possible association with the ra- 
dio source MGPS-2 Jl 10149-610104, which has a flux of 
24.2 ± 4.8 mJy at 843 MHz. They suggest that IGR J11014- 
6 1 03 may be a pulsar wind nebula (PWN) produced by a high- 
velocity pulsar but conclude that X-ray observations with bet- 
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ter angular resolution are required to establish the nature of 
the source. 

Here, we report results of an observation with the Chan- 
dra X-ray Observatory, which provides the improved angular 
resolution. 

2. CHANDRA OBSERVATION AND DATA REDUCTION 

A Chandra observation of IGR Jl 1014-6103 occurred on 
2011 September 6 (ObsID 12420). We used the Advan ced 
CCD Imaging Spectrometer (ACIS. iGarmire et al.l [200 3) as 
the detector with the INTEGRAL position of IGR J11014- 
6103 at the ACIS-I aimpoint. The observation started at 
1.43 hr UT, and a net exposure time of 4,991 s was obtained. 

The data were initially processed at the Chandra X-ray 
Center (CXC) with ASCDS Version 8.4. After obtaining 
the data from the CXC, we performed all subsequent pro- 
cessing with the Chandra Interactive Analysis of Observa- 
tions (CIAO) version 4.3.1 software and Calibration Data 
Base (CALDB) version 4.4.6. We used the CIAO program 
chandra_repro to produce the "level 2" event lists, which 
we analyzed as described below. 

3. ANALYSIS AND RESULTS 

An inspection of the 0.3-lOkeV ACIS image indicates the 
presence of the complex source that was previously reported, 
and the point source (PS) and North-East (NE) extension are, 
by far, the brightest sources in the field-of-view. Figure [T] 
shows NE in far more detail than in previous observations, 
indicating that the emission may, in fact, have a conical or 
cometary shape with its narrow end close to PS. 

3.1. The Point Source 

We determined the Chandra position for PS using the CIAO 
routine wavdetect, and it is R.A. = 1 1''01'"44".96, Decl. = 
-61°01'39".6 (equinox 2000.0) with a 90% confidence un- 
certainty of 0".64. Although the previously reported XMM- 
Newton position is marginally compatible with the ~'0".l 
(1-cr) near-IR position of 2MASS Jl 1014532-6101383, the 
Chandra position is 3".0 away, which rules out the associ- 
ation. Furthermore, we searched all of the source catalogs 
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Fig. 1.— Chcmdm/AClS 0.3-10 keV image of IGR Jl 1014-6103. There is 
a point source (PS) and an extension to the North-East. The image has 0".98 
pixels and has been smoothed with a Gaussian kernel with a radius of 2". 95. 

available in the VizieR databas43 with a search radius of 2".0, 
and the XMM-Newton source is the only one that appears. Fi- 
nally, we obtained near-IR images from 2MASS and optical 
images from the Digitized Sky Survey. Although there is no 
evidence for near-IR or optical sources at the Chandra po- 
sition of PS, the nearby 2MASS source is relatively bright 
(J = 11 .2, / = 13.0), making it difficult to obtain precise upper 
limits on the near-IR and optical brightnesses without higher 
quality images. 

We also produced a radial profile for PS to determine if it is, 
in fact, a point source, and the profile is shown in Figure|2] We 
used Chandra Ray-Tracing software (ChaRT) and the Chan- 
dra simulator ("MARX") to calculate the point spread func- 
tion (PSF) at the position of PS. The inputs to ChaRT were 
the off-axis angle (3'.64), the azimuth (330.82°), and the best 
fit spectrum for PS described below. Figure |2] shows that PS 
is consistent with being a point source. In fact, the simulated 
profile is slightly higher than the measurements for PS, and 
this is likely due to inaccurac ies within ChaR T or MARX as 
has also been found by Bogdan ov et alj ( 1201 lb . 

We extracted an energy spectrum for PS using the CIAO 
program specextract. We used an aperture with a ra- 
dius of 5" for the source spectrum, and a nearby rectangu- 
lar source-free region for the background spectrum. We ex- 
tracted 126 counts in the source region in the 0.3-10 keV 
band, and the background estimate (scaled to the size of 
the source region) is 0.8 counts. We fitted the energy spec- 
trum using the XSPEC software package. A model con- 
sisting of a power-law with interstellar absorption provides 
a good description of the spectrum. The power-law pho- 
ton index is F = 1.5;I;o 5 (90% confidence uncertainties are 
given here and below), and the unabsorbed 0.3-10 keV flux is 
(8.6!^ ^) X 10"'^ ergcm-2 s"'. UsinglWilms. A llen & McCrav 
(2000) abundances and iBalucinska-Church & McCammon 
(|l992) cross sections, we obtained a column density of A^h = 
(l-fi^^ i) X 10^^ cm"-. This is cons istent with the pr evious 
value measured with XMM-Newton (iPavan et alJl201 II) . 

We also checked for long and short-term variability of 

^ See |http://webviz.u-strasbg.fr/viz-binA^izieR-4| 
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Fig. 2. — The radial profile of PS (solid line). The dotted line is the Chan- 
dra point spread function from "ChaRT" at the position of PS and with an 
energy spectrum like that of PS. 

PS. Using the spectral shape described above, the Chan- 
dra observation indicates an absorbed 2-lOkeV flux of 
5.5 X 10~'^ergcm~^ s"', and the flux obtained in the same 
band using the XMM-Newton observation is (6.2;|;2 5) x 
10~'^ergcm"^ s"' (iPavan et al.ll201 Ih . Thus, there is no ev- 
idence for flux variations bet ween the 2003 and 2011 obser- 
vations. Using XMM-Newton, IPavan et all (1201 Ih find no in- 
dication of variability on time scales of seconds to hours. We 
searched for variability with Chand ra down to a time scale of 
50 s using the CIAO tool glvary (iGregorv & LoredollI992l) 
and did not detect variations. 

3.2. The North-East Extension 

We extracted an energy spectrum for NE from a 24 "-by- 
46" rectangular region rotated to a position angle of 45°. We 
chose the region to include all of the counts from NE and to 
exclude the counts from PS. We extracted 261 counts in the 
source region in the 0.3-lOkeV band, and the background es- 
timate (scaled to the size of the source region) is 1 1.7 counts. 
A model consisting of a power-law with interstellar absorp- 
tion provides a good description of the spectrum. The power- 
law photon index is F = 1.8 ±0.4 and the column density is 
A^H = (0.8!*;:^) X 1022cni-2. The unabsorbed 0.3-lOkeV flux 
is (1.3;!;o2) ^ 10~'^ergcm"^ s~'. The spectra l param eters are 
consistent with those found by IPavan et al] ( 1201 II) . but the 
flux measured by Chandra (the absorbed 2-IOkeV flux is 
6.4 X 10"'^ ergcm"^ s~') is a factor of ^1.7 higher It is pos- 
sible that this difference reflects a true change in flux, but it 
may also be related to the fact that PS and NE are not fully 
resolved by XMM-Newton. 

Profiles of detected counts for PS and NE are shown in Fig- 
ure |3^. We extracted counts from a 24"-by-80" rectangular 
region rotated at a position angle of 45° that includes both PS 
and NE. We produced profiles of the counts along the long 
axis of the rectangle in two energy bands, 0.3-lOkeV and 2- 
lOkeV. While the number of 2-lOkeV counts is close to half 
the number in the full energy band, the profiles are not signif- 
icantly different. The profile for NE peaks between 15" and 
30" from PS. 
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Fig. 3. — (a) Counts profiles for PS and NE in two energy bands: the full 
0.3-10 keV band (black) and the 2-10 keV band (grey). The horizontal axis 
runs through PS and NE at a position angle of 45°, and the counts in the 
perpendicular direction within ±12" of the axis are summed, (h) Unbinned 
0.3-10 keV image where each dot corresponds to the location of an X-ray de- 
tection. The dotted, solid, and dashed lines represent the predicted locations 
of the bow shock for 5so = 0".4, 1 " , and 2" , respectively, assuming i = 0° . 

3.3. The X-ray Tail 

In addition to PS and NE, the X-ray "tail" that was 
previously reporte d from the XMM-Newton observations in 
iPavan et al] (1201 Ih is also detected (see Figure|4]l. The length 
of the tail is several arcminutes, and it is nearly perpendicular 
to the NE extension. 

4. RADIO PULSAR SEARCH 

We searched for pulsations from IGR J 1 1 1 4-6 1 03 on 20 1 1 
July 27 using the CSIRO Parkes radio telescope. We ob- 
tained a 5.6 hr observation with the 1.4 GHz multibeam re- 
ceiver and sampled every 1 ms each of the 96 frequency chan- 
nels spanning 288 M Hz of bandwidth (for f urther data ac- 
quisition details, see [Manchester et"al]|2001l) . We analyzed 
the data following standard pulsar search techniques using 
PRESTO (Ransom 2001), but did not identify a convinc- 
ing pulsar candidate. We dedispersed the data up to DM = 
1280pccm"^, which is twice the maxim um predicted in this 
direction by the ICordes & Laziol (12002b model for the Galac- 
tic distribution of free electrons. The sensitivity limit of this 
observation, assuming a pulsar duty cycle of 10%, is S\a = 
0.05 mJy, corresponding to a luminosity (Li 4 = S\Ad^) limit 
of <0.05 mJy kpc^ at a distance of 1 kpc and <5 mJy kpc^ at 
d = 10 kpc. If the source is at the larger distance, the non- 
detection is not surprising since approximately half of all 
young pulsars detected in the radio have luminosities below 
SmJykpc^ (see lCamilo et allbOOQl) . 

5. DISCUSSION 

The results presented above are consistent with (but do not 
prove) the hypothesis that PS is a pulsar. The Chandra posi- 



tion shows that the source does not have a bright optical/IR 
counterpart, which is consistent with a pulsar since isolated 
neutron stars are very faint in the optical and IR. The stabil- 
ity of the long-term flux and the power-law spectrum are also 
expected if the source is a pulsar. 

The Chandra observation provides new information about 
the NE extension and its relationship to PS. In the scenario 
where PS is a pulsar, th e natural explanatio n for NE is that it is 
a PWN as suggested bv lPavan et al.l (1201 Ih . For some PWNe, 
a bow shock is created due to the motion of the pulsar through 
the ISM. In the case of IGR Jl 1014-6103, the morphology of 
NE suggests that the pulsar may be moving in the South- West 
direction at high speed. Figure |4] provides evidence in favor 
of this interpretation as it shows that NE is elongated in a 
direction that is co-aligned with PS and the center of the SNR 
MSH 11-61 A. Thus, a strong possibility is that the compact 
object that is powering IGR Jl 1014-6103 was produced when 
the supernova occurred. 

MSH 1 1-61 A has a centrally bright X-ray morphology, but 
the location of the compact object produced in the supernova 
is unknown. The X-ray emission from the SNR is ther mal 
rather than being powered by a PWN ( Slane et alll200 2V It 
has been suggested that PSR J 1 105-6107, which is 22' from 
the c enter of the SN R, may be associated with MSH 11- 
61A (i Kaspi et aT][T997.). Applyin g SNR evolution models to 
MSH ll-61A. ISlaneetalJ dSool argue that the SNR is (1- 
2)x 10'*yr old and is at a distance of 8-1 1 kpc. The thermal 
conduction and cloudy ISM models would imply transverse 
velocities f or PSR Jl 105-610 7 of 4,500 and 5,300 kms"', re- 
spectively (iSlane et al.ll2002h . In contrast, the known distri- 
bution of pulsar velocities includes groups near 175 km s"' 
and 700km s-i ( ICordes & Chernoffi >1998). The fastest pulsar 
with parallax and pr oper motion measurem ents has a veloc- 
ity of 1,100km s"' (IChatterjee et all 120051) . Although there 
is evidence for pulsar or compact object velocities in ex- 
cess of l,100kms-i (Hui& Becker 20061: IWinkler & Petrd 
^2mh , > 1,000 km s"' ^ getal. 2011), and > 1,400 km s"' 
(.Lovchinskv et al.ll201 Ih . the velocities implied by an associ- 
ation between PSR Jl 105-6107 and MSH 1 1-61 A are likely 
too high. 

IGR Jl 1014-6103 is 11'.9 from the center of MSH 11-61 A, 
and scaling from the calculations previously carried out for 
PSR J 1 105-6107, we determine transverse velocities of 2,400 
and 2,900km s"' for the thermal conduction and cloudy ISM 
models, respectively, if IGR Jl 1014-6103 was formed in the 
supernova event that produced MSH 1 1-61 A. Among known 
compact objects associated with SNRs, these values are 
only rivaled by the recent est imate of 1,400-2,600 km s"' for 
XMMU J172054.5-372652 (iLovchinsky et alJiMlh . Thus, 
these two cases indicate that supernova models that can pro- 
duce >2,000km s"' kick velocities may be required. 

5.1. The PWN Bow Shock Interpretation and Implications 

If the interpretation that IGR Jl 1014-6103 is a PWN and 
that NE is emission from a bow shock is correct, then the fact 
that PS is a point source along with the morphology of NE 
provide a constraint on the angular distance between the pul- 
sar and the standoff shock, S^o = Rso cos // d, where /?so is the 
standoff shock radius, / is the inclination angle of the pulsar 
velocity with respect to the plane of the sky, and d is the dis- 
tance to the source. The standoff shock location is set by the 
ram pressure balance between the pulsar wind and the pulsar 
passing through the interstellar medium (ISM). The shape of 
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Fig. 4. — Composite image of tlie supernova remnant MSH 11-61 A near IGR Jl 1014-6103. The large image is from XMM-Newton (0.5-10 keV), while the 
inset is from Chandm-AClS (0.3-10 keV): both were rebinned and smoothed. The blue circle is the 4'. 3 INTEGRAL error circle. It is shown via the red arrow 
that sources PS, NE, and the center of the SNR, marked with a cyan square, are co-aligned. The 4' X-ray tail to the North- West of PS can also be seen. 



the entire bow shock in space is determined by R^o according 
to 



for the standoff shock radius becomes 



^(</>)=^l/3(l- 



smr 



tan( 
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where (p is the angle with respect to the direction of motion, 
and R is the distance to the shock along each angle ( Wilklfl 
|2003). While the observed shape of the bow shock will de- 
pend on /, here, we make a rough estimate for S^o using the 
/ = 0° case. Figure [3J5 compares the calculated shape of the 
bow shock to the Chandra data for three values of 6so- The 
data are best described by a value close to 1 ", and the lower 
and upper li mits are approximatel y 0".4 and 2", respectively. 
Following lCaraveo et al.l d2003h , we derive 



Sso = 0".266 



cos^ / 



Vf.3 diQ no. 



(2) 



where v,_3 is the transverse velocity of the pulsar in units of 
10^ km s"', dw is the distance to the pulsar in units of lOkpc, 
£36 is the rotational energy loss of the pulsar in units of 
10^^ ergs"', and no.i is the particle density of the ISM near 
the source in units of 0. 1 cm"^. 

From the spectral fits described above, the 0.3-10 keV 
unabsorbed flux of PS and NE combined is 2.2 x 
10"'^ erg cm"^ s"', which corresponds to a luminosity of 2.6 x 
lO-'"^ c/igergs"', and as rj is the radiative efficiency of the 
PWN, £36 = 0.026t/^Q77~'. We obtain another constraint by 
considering that if NE is a PWN, then the emission has a syn- 
chrotron origin and will decay as the electrons age on a time 
scale of fs_3, where this is the exponential e-folding time for 
the X-ray flux in units of 10^ yr Thus, the ~1' transverse 
length of NE impHes that V;^ = 2.8t^iof~3, expression 



6so = 0".015 



cos ; f5_3 
dw njf 771/2 



(3) 



A significant caveat is that the synchrotron lifetime con- 
straint assumes that any velocity imparted to electrons 
in the PWN by the pulsar is small relative to the pul- 
sar space velocity. If NE has a ma gnetotail component 
(iRo manova, Chulskv & Lovelacell20 05'), then the ejection ve- 
lociti es could be high enough to invalidate this assumption 
(Kai-g dtsevetaLll2008l) . 

Figure |5] shows the distance-^so parameter space for the 
case of / = 0°, f.5,3 = 1, and no.i = 1 for several differ- 
ent values of 77 that span values measured for '^60 PWN 
^Kargaltsev & Pavlov 2010). The value of 3 = 1 is obtained 
in the lCaraveo et al.l (i2003) study, an d it is what would be pre - 
dicted for a ~10/LtG magnetic field jKargaltsev et alj|2008l) . 
The value of no 1 = 1 is consistent with the intercloud ISM 
density foun d in the vicin ity of MSH 1 1-6 1 A using the cloudy 
ISM model (ISlane et al.l l2002). The Chandra constraints on 
Sso suggest a value of 77 10"^ if the distance is ^2kpc and 
values of 77 ^ lO'^-lO"'* if the distance is ^lOkpc. Although 
we cannot use this information to formally constrain the dis- 
tance, the larger distance is slightly favored since PWNe with 
rj = 10"^- 10"^ are mu ch more common than PWN e with larger 
radiative efficiencies dKargaltsev & Pavlovll2010l) . 

5.2. X-ray Tails from PWNe 

Whether IGR Jl 1014-6103 is related to MSH 11-61 A 
or not, if PS is a pulsar moving in a South-West direc- 
tion as implied by NE being a bow shock, then the X-ray 
tail extends ~4' in a direction that is ^perpendicular to 
the motion of the pulsar. A chance alignment seems un- 
likely, especially since similar X-ray tails have been seen 
from other PWNe. Examples include the Guitar Nebula 
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Fig. 5. — Parameter space for the distance to IGR Jl 1014-6103 and the 
standolf shock radius ((5so)> if NE is indeed a bow shock. From the Chandra 
image, we find that values of 5so between 0".4 and 2".0 are allowed (for 
i = 0°). The solid lines show the calculated values for a range of pulsar 
efficiencies (rj) for the case where i = Q° ,no \ = 1, and 3 = 1. 

(PSR B2224+6 5) and PSR Jf0357+3205. which have 2 ' and 
9' X-ray tails ( Hui & BeckelllOOTHDeLuca et al.ll20TTl) . The 
Guitar Nebula is noteworthy because its X-ray tail points in a 
direction that is 118° from the measured proper motion direc- 
tion for this pulsar (Hui & Becker 2007; Hui et al. 201 1). 

While these examples support the PWN nature of 
IGR Jl 1014-6103 as discussed by iPavanetal] ( 1201 Ih and 
also the possibihty that the IGR Jl 1014-6103 X-ray tail may 
extend in a direction perpendicular to the pulsar's direction of 
motion, it is worth pointing out that the size of the X-ray tail 
may be an argument against a ^10 kpc distance. While the X- 
ray tails for the other pulsars have lengths near 1 pc, a 10 kpc 
distance for IGR Jl 1014-6103 would imply a --10 pc length 
for the X-ray tail from this source. 



6. SUMMARY AND CONCLUSIONS 

This work provides further evidence that IGR Jl 1014-6103 
is a PWN. The Chandra position for PS eliminates the previ- 
ously suggested 2MASS counterpart, and the available near- 
IR and optical images do not indicate the presence of another 
counterpart, which is consistent with the faint levels of opti- 
cal/lR emission expected from a pulsar Most significantly, 
the Chandra observation has provided new information about 
the morphology of the extended source to the North-East of 
PS (called NE in this work), and we argue that the most likely 
interpretation is that NE is shaped by a bow shock caused by 
PS moving in the South- West direction at a high rate of speed. 

Given the alignment of PS, NE, and the supernova remnant 
MSH 11-61 A, we suggest the possibility that IGR J11014- 
6103 is a pulsar that was formed by the supernova event that 
produced MSH 11-61 A. If correct, this would imply that 
IGR Jl 1014-6103 is at the same distance as the SNR, 8- 
1 1 kpc, and also that the supernova imparted a kick velocity of 
at least 2,400-2,900kms"'. While we show that the Chandra 
constraint on the standoff shock makes a scenario with a large 
distance and a high velocity plausible, caution is required due 
to uncertainties about the nature of the source, its orientation 
in space (/), and its environment. One possible counterargu- 
ment to IGR Jl 1014-6103 having a distance compatible with 
the SNR is that this would imply a lOpc-long X-ray tail. 

We acknowledge useful discussions with J. Halpem, 
V. Kaspi, and S. Boggs. Support for this work was provided 
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of the Australia Telescope, which is funded by the Common- 
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